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Abstract 

The large elliptic flow observed at RHIC is considered to be evidence for almost perfect liquid 
behavior of the strongly coupled quark-gluon plasma produced in the collisions. In these pro- 
ceedings we present a two parameter fit for the centrality dependence of the elliptic flow V2 scaled 
by the spatial eccentricity e. We show by comparing to viscous hydrodynamical calculations that 
these two parameters are in good approximation proportional to the shear viscosity over entropy 
ratio rjl s and the ideal hydro limit of the ratio vi/e. 



1. Introduction 

The goal of the ultra-relativistic nuclear collision program is the creation and study of a new 
state of matter, the quark-gluon plasma. The azimuthal anisotropy of the transverse momentum 
distribution in non-central heavy-ion collisions is thought to be sensitive to the properties of this 
state of matter. The second Fourier coefficient of this anisotropy, V2, is called elliptic flow. For a 
recent review see ||T|. 

In ideal hydrodynamics V2 is proportional to the spatial eccentricity with a magnitude which 
depends on the Equation of State EoS. This spatial eccentricity is defined by 

„ ^ - x^) 

where x and y are the spatial coordinates of the colliding nucleons in the plane perpendicular 
to the collision axis and where the brackets denote an average. In practice e is not a measured 
quantity but obtained from model calculations, using Glauber or Color Glass Condensate (CGC) 
models, for instance. 

The ratio V2/e versus particle density is a sensitive gauge to test if the system approaches 
ideal hydrodynamic behavior [2]. It was observed that this ratio reaches the expected ideal 
hydrodynamic values only for the more central collisions at the highest RHIC center of mass 
energy [3 , 4J which indicates that certainly for non-central collisions, as well as at lower energies, 
and away from mid-rapidity the elliptic flow contains significant non-ideal hydro contributions. 

Much of this discrepancy can be explained by incorporating the viscous contribution from the 
hadronic phase |5 6, 7j. However, we expect that also the hot and dense phase must deviate from 
an ideal hydrodynamic description. Kovtun, Son and Starinets (KSS) lO, showed that conformal 
field theories with gravity duals have a ratio of shear viscosity rf to entropy density s of, in 
natural units, tjI s - IjAn. They conjectured that this value is a lower bound for any relativistic 
thermal field theory. In addition, Teaney Q pointed out that very small shear viscosities, of the 
magnitude of the bound, would already lead to a significant reduction in the predicted elliptic 
flow. 
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Based on the centrality dependence of V2/e, the magnitude of 77/s for the created system has 
been estimated recently from a transport theory motivated calculation 1.10. ,1 IJ and from viscous 
hydrodynamical calculations lfT2lfT3l . Both approaches have their merits and drawbacks. 

In these proceedings we explore how well a parameterization can be used to estimate 77/5 as 
well as the ideal hydrodynamical limit of vi/e which is closely related to the EoS. 



2. Simple Parameterization 

We use the parameterization from |l2l[T0l which is defined by 



where S is the transverse area of the collision region and h and B are the two free parameters in 
the fit. The parameter h corresponds to the ideal hydro limit of V2/e and B is proportional torj/s. 

Figure [T] shows how the parameterization behaves for two diff'erent values of the ideal hydro 
limit (the dashed line represents the harder EoS) and two different values of ri/s (the full line 
represents the smaller rj/s). The effect of the EoS is clearly seen in the magnitude of V2/e in 
Fig. [T] and the value of rj/s is reffected by the change in this magnitude versus l/SdN/dy (for 
rj/s - the magnitude will be constant). The magnitude of rj/s is easier to quantify if one 
plots vijhe, as is done in Fig. |2] A larger deviation from unity at fixed value of 1 /SdN/dy then 
indicates a larger 77/ s. 

To test if this simple parameterization does describe a state of the art viscous hydrodynamical 
calculation we fit the calculations from Luzum and Romatschke lfT4ll . Figure |3] shows that Eq.[T| 
well describes results from viscous hydrodynamical calculations, done with three different values 
of rj/s and two different parameterizations of the spatial eccentricity (Glauber and CGC). As 
expected, V2 is to good approximation proportional to the initial spatial eccentricity. In addition, 
it is seen that the deviation of V2/e from unity at a given 1 /SdN/dy increases for larger values of 
77/s. 

Figure [4] shows V2/hs from viscous hydrodynamical calculations lfT2l [T3l [141 done by differ- 
ent groups using the same set of values of rj/s but different parameterization of the EoS and s. 
The value of e is that used in the hydrodynamical calculations while the value of h is obtained 
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Figure 1: The dependence of vi/e versus transverse density Figure 2: The dependence of V2/hs versus transverse den- 
of equation^for two values of h and two values of i. sity of equation^for the same parameters as Fig.^ 

2 




20 40 60 80 100 120 "o 20 40 60 80 100 120 

1/S dN/dy 1/S dN/dy 



Figure 3: A fit of viscous hydrodynamical model results Figure 4: Comparing viscous hydrodynamical calculations 
using CGC and Glauber initial eccentricities with Eq.[T] of different groups with the fit 

from the fit. We conclude that our parameterization yields curves that depend on the value of 
rjl s but are roughly independent of the EoS and e. However it turns out that if the EoS is very 
different (e.g. not incorporating a phase transition) this scaling does break down (not shown). 

Using Eq. [T] we can now compare the various viscous hydrodynamical results with data and 
estimate the value of rj/s. Since the value of s is not known we take the eccentricity calculated 
assuming CGC IfTSi or Glauber (wounded nucleon) initial conditions as two extremes. It is seen 
from Fig.|5]that, assuming the CGC initial conditions, the STAR data is well described with twice 
the KSS bound, tjI s - 21 An. Using the Glauber initial conditions, however, the STAR data is 
not described within the range of rjl s currently used by the viscous hydrodynamic calculations. 
From the deviation from unity one can estimate that the corresponding value of rjl s would be 
approximately four times the KSS bound. Using the CGC or Glauber initial conditions we find 
for the ideal hydro limit of V2/e the value 0.20 ± 0.01 and 0.36 + 0.07, respectively. 

For the CGC initial conditions the value of h approximately matches the EoS used by Luzum 
and Romatschke |13|). This is illustrated in Fig. |6] where the centrality dependence of vo [14J is 
well described by CGC initial conditions, a value oih ^ 0.2 and 77/s = 2/47!-. 

Using viscous hydrodynamics with these CGC initial conditions, EoS, and magnitude of 
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Figure 5 : Comparing viscous hydrodynamical calculations Figure 6: A direct comparison of viscous hydro calcula- 
with STAR data. tions with PHOBOS data (from fT4l l 
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Figure 7: V2 from STAR (approximately corrected for nonflow) compared Figure 8: V2 from STAR as function of 
to viscous hydrodynamical calculations (from U3J ). transverse momentum and centrality. 

rj/s, the transverse momentum dependence of V2 is also well described, as shown in Fig. |7] 
The figure illustrates that the pt dependence is very sensitive to the viscous correction such that 
larger corrections decrease the magnitude of v'2 and shift its maximum to lower p^. Figure [8] 
shows the centrality dependence of V2{pt) where one clearly observes that the deviation with 
rj/s = increases from central to peripheral collisions and that the peak position shifts to lower 
Pi, consistent with larger viscous effects. 

3. Conclusions 

We have shown that a simple parameterization can describe the centrality dependence of 
va/e. When compared to viscous hydrodynamical calculations such a parameterization yields an 
estimate of rj/s. We find that the current RHIC data is described well by a spatial eccentricity 
based on CGC initial conditions, a soft EoS with va/e ~ 0.2 and t]/s twice the KSS bound. 
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